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Abstract 

The large N phase transition point is investigated in the heat kernel on 
the U(N) group with respect to arbitrary boundary conditions. A simple 
functional relation is found relating the density of eigenvalues of the bound- 
ary field to the saddle point shape of the typical Young tableaux in the large 
./V limit of the character expansion of the heat kernel. Both strong coupling 
and weak coupling phases are investigated for some particular cases of the 
boundary holonomy. 
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1 Introduction 



Heat kernels on classical groups are widely used in mathematical physics, especially 
in quantum field theory [p]] [[|. For example it serves as a building block (disc parti- 
tion function with an arbitrary boudary holonomy) for calculations of the partition 
function of Yang-Mills theory on 2D manifolds with various topologies @ I! • 
The large N limit of the heat kernel is of special interest since it is related to the so 
called planar limit in some quantum field theories, QCD 2 being the most interesting 
among them |||||. This limit not only can simplify the analytical investi- 



gation of a theory but can also give rise to new qualitative phenomena, like large TV 
phase transitions [TT[ . 



Our task in this paper is to calculate, as explicitly as possible, the heat kernel for 
the U(N) group (as an example) in the large N limit. The method will be closely 
related to the recent calculation of the partition function of Yang Mills theory on the 



two dimensional sphere [11 



To give a physical taste to the heat kernel let us derive the heat kernel as the 
partition function of a simple quantum mechanical model; the one dimensional prin- 
cipal chiral field with initial and final holonomies specified . The functional integral 
representation of the partition function is 

Z(T, 9 , Br) = J [Vg(t)) H e mr Jo d ^- ls t9, (l) 

where [Dg(t)]H is the Haar measure on the unitary group U(N), and the holonomies 
at the end points of the time interval are g(0) = e ld ° and g(T) = e ldT (with e te ° and 
e ldT the eigenvalues of g(0) and g(T)). Due to the invariance of the Haar measure at 
each moment of time, we have, after the introduction of a new variable (connection) 
A(t) = ig~ x d t g{t) the following representation of Z(T,9 ,9t)' 

Z(T } 6 0} 9t) = J VA(t)e- Ntr Jo MA V 2 S^Te'Io dtA ^]^, e WT ). (2) 

Using the character expansion of the group ^-function: 8(U, U') = J2r Xr(U)xr{U')*, 



and the fact that the character is the trace of the matrix element in a given representa- 
tion R: xn(Texp(i J Adt)) = tr R [T exp(i J A a r^dt)] (rg is the ath generator of U(N) 
in the Rth irreducible representation), we arrive after a simple gaussian integration 
in A(t) (independant at each moment t) at the conclusion that 

Z(T,e ,6 T ) = J2^ TC2{R) XR(e td0 )XR(e* dT ) (3) 

R 

where C 2 (R) is the second Casimir of the group U(N) in the representation R. This 
is also equivalent to the partition function of two dimensional QCD on the cylinder 
with the boundary holonomies specified at either end of the cylinder by the two 
distributions e l9 ° and c i6t . For the rest of this article we restrict our attention to the 
case where one of the holonomies is the unit matrix (two dimensional QCD on the 
disc), and the heat kernel takes the form 

Z(T,e) = Y,dne- TC2{R) XR(e i6 ). (4) 

R 

In section(2) we derive an explicit formula for the partition function for the heat 
kernel with fixed boundary field. The solution is written directly in terms of the 
eigenvalues of the boundary field. The calculation is valid only in the weak coupling 
phase. In section(3) we find a simple inversion relation between the density of eigen- 
values of the boundary field and the density of heighest weights labeling the saddle 
point Young tableau. In some important particular cases of boundary holonomy, 
such as the semicircle distribution, it is shown to be explicitly solvable. From this it 
is then trivial to calculate the strong coupling transition point. We then discuss in 
section(4) how to calculate the partition function in the strong coupling phase. An 
alternative approach is discussed in the appendix. The appendix also contains an 
explicit solution for the heat kernel action illustrating the procedure developped in 
section(3) for a nontrivial boundary field (the density of the #'s being given by an 
inverted semicircle distribution a{9) = — , 1 ). 
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2 Weak coupling phase 



To study the action in the large N limit we proceed as in JTT|. The sum over repre- 
sentations R for the group U(N) is the sum over all Young tableaux labeled by the 
components of heighest weight n\, tin- The rij are integers, and the sum over 

representations is given by summing over all {rii} satisfying the inequality: 

oo > rii > n 2 > .... > n N > — oo. (5) 

With this labeling for the representation, the explicit expressions for the dimension, 
second Casimir and character are: 

d R = U i>j (l-^ i ), (6) 
C 2 (R) = Eli ni(m -2i + N + l), (7) 
Xn(e* s ) = kA A{ l„ } \ (8) 

where the determinant is of the matrix whose (k, I) element is e l ( n k-k+N)9i anc j ^.(e 16 ) 
is the Vandermondt determinant A (a;) = n^i^i ~ x j)- 

It is convenient to introduce a new set of labels, {h^} defined by 

h i = -n i + i-{N+ 1)/2, (9) 

which satisfy the inequality 

h i+1 > hi. (10) 

The heat kernel is then given by 

Z= &(h)e-^?=i h Uet kl (e- ihk9l )f(e) (11) 

{hi} c onstrained 



where 



f(9) = ^e-^-^^\ (12) 



A(e 



and the sum over the {hi} is an ordered sum subject to the inequality (pT0|) . The 
determinant consists of iV! terms of alternating sign. By interchanging the {hi}, 



each of these terms can be made equal to e^ 1 ^^ hi6 '^, the ordered sum becomes 
unconstrained, (each of the N\ terms contributing a different segment of the sum 
over all {hi}) and the partition function becomes: 

Z=J2 A ( h ) e-Efei^-KMO (13) 

{hi} 

In the large N limit we redefine the {hi}, dividing them by N, and then treat them 
as continuum variables so that the multiple sum becomes a multiple integral: 

oo 

Z = J (j[dh t ) A(h) e-NEi^hiOi (14) 

— oo * 

It is important to realise that by changing the sum into an integral we may have lost 
the inequality (fLOf). The expression for Z we now obtain is therefore only valid in the 
weak coupling phase where (^) is satisfied. By shifting the {hi} we can complete 
the square and do the gaussian integral: 

Z= /(ilidfc) A(h-^9)e- N ^ Th i + ^i f(6) (15) 
A(i0)e-&£i«? T- N2 ' 2 f(6). (16) 

The last step is accomplished by noting that the answer is a completely antisymmetric 
polynomial of the 9i of the same degree as the Vandermondt determinant, ie., it is 
the Vandermondt determinant up to a constant. We have thus restored the partition 
function in the weak coupling phase ||. We will see in the next section that the weak 
coupling solution will be valid only for small enough T, due to the constraint flT0|). 

3 Strong Coupling Transition Point 

It is expected that in the large N limit the integral above is dominated by the 
contribution from a single set of {hi} ie., a single representation. In terms of the 
continuum variables, this corresponds to a density p(h) = jf§K f° r the distribution of 
Tt's. As is discussed in [fLl|| , for this distribution to correspond to an allowed Young 
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tableau p(h) must satisfy the constraint p(h) < 1, following directly from (|H]). The 
point at which p(h) = 1, for some h, corresponds to the strong coupling transition 
point beyond which the weak coupling solution is no longer valid, and the strong 
coupling methods discussed in the following section have to be used. The rest of this 
section is entirely devoted to finding the distribution p(h) corresponding to the weak 
coupling solution, to allow us to test for this condition. 

To find the density we calculate the resolvent, H(h), which is defined to be: 

ff w = G£rhjl- (17) 

where h lies in the complex plane and the average is taken with respect to h'. In the 
large N continuum limit this is given as the integral along the support of p(h) 



H(h) = J dh- 



h-h' 1 

The discontinuity across the cut of the resolvent gives 27ri times the density, p(h). 
To proceed we take the Fourier transform of each term in the sum (|TT|): 



N ■ 00 



1 i r 

H ( h ) = Jf E g J dk ^me- lkh (e^ k ) hl (19) 

The expectation value, (e ih '^ k )h', can be calculated from the formulae (|14D and (|16|) 
for Z by realising that the term e lhlk can be absorbed into the definition of the 9. 
We therefore define a new distribution of #'s which is identical to the old distribution 
except for the zth 9 which is shifted by —k/N; 

Oj = e 3 - pij. (20) 

This allows us immediately to write down the result 

_ JV_ g2 i 

n = = n (i - m^y^- 2 ^ (2D 

Substituting this into equation (|19D we notice that the sum can be replaced by a 
contour integration in 9 (encircling all the 9 eigenvalues) if we also unrestrict the 



product, so that the product is over all j: 

H{h) _ j- dk ^: hj ^ _ _t_ K *«^) m 

We can then exponentiate the unrestricted product and keep only the lowest order 
in 1/N to arrive at 

1 I jn f°° ji. AL\ -cosk(h + i^-iQ(9)) 



H(h) = — i d6 / dk (e~oT) v , 2 (23) 

2n J Jo k 

where 

is the resolvent for the distribution of 9. The apparent divergence of the k integral 
about k = is canceled by the contour integration in 9 (we can add any non 9 
dependant function of k into the k integrand, eg. (e"^ 7 ) 00 ^ , since it will be 
ignored by the contour integral). The gaussian exponential term suppresses the 
divergent behaviour of the cos (it's argument is complex) and permits us to do the 
integral unambiguously in the large N limit, the result being: 

H{h) = ^f c d9 \og[h + i A - iG(9)\ (25) 

where the contour C goes anticlockwise around the cut of the resolvent 0(9). We 
would like to thank Matthias Staudacher for simplifying this derivation. 

In the appendix we show how to find an analagous formula (the roles of H(h) and 
©(#) are interchanged) using either Migdal's method [HJ or some results developped 
by Matytsin |4|. 



It is more convenient to study the derivative of H(h): 

H'(h) = — I d9- g- — (26) 

In this case the cut point (s) of the logarithm becomes a pole(s) and the contour can 
be inflated out from the cut to catch both the pole(s) and also a contribution from 
infinity: 

H ' {h) = 2T -L[^-mem (27) 
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where the 6(h) satisfy 



h + i e M-ie(e(h)) = o (2? 



and the sum over the zeros means the sum over all the 6(h) which satisfy equation(^) 
Taking the derivative with respect to h of equation(|28|) gives an identity which allows 
us to reduce equation (|27|) to 



H'(h)=2T + iJ2 8'(h) (29) 

zeros 

Integrating up with respect to h we arrive at the very simple formula 

H(h) =2Th + iJ2 6(h) (30) 

zeros 

Some general comments are now in order. Firstly, we only need to solve the inversion 



equation(p8|) in some neighbourhood of h, as analytic continuation then defines H(h) 
over the whole complex plane. The cutpoints are completely specified, although the 
actual path of the cuts between the cutpoints maybe ambiguous, in which case they 
are chosen to lie on the real axis. For densities a(6) which are non-singular, ie. finite 
on all of their support, we can choose the neighbourhood of h to be far from the 
origin so that there is only one solution of equation (|28|) . In this case we have 

H(h) = 2Th + i6(h) = 2TiQ(6(h)) (31) 



where the last equality follows directly from (^Sj). Studying the large h behaviour of 
the right hand side we find that H(h) behaves as jr for large h as is required. For this 
reason we have neglected the possible constant term in equation(|30|) coming from the 
integration. 

We are now in a position to find the strong coupling transition point, ie. to test the 
inequality p(h) < 1. If the distribution of Tt's reaches it's maximum at h = we can 
derive a very simple expression for the strong coupling transition point. Setting h = 
in equation (|3~T1) we have at the strong coupling transition point Q(6(0))) = ±7r/2T, 



which implies (from equation(|28|)) that #(0) = ±7r which leads in turn to the simple 



expression for the strong coupling transition point: 



T < = 4o (32) 



We would like to thank Andrei Matytsin for this comment [12 . 



3.1 An example: The semicircle distribution 

To illustrate the method discussed above we present the simplest example, the case 
where the spectrum of 0's is given by a semicircle distribution: 



a(6) = —V^¥ and 6(0) = -(0 - Vff^c) (33) 

7TC C 



Inverting equation (|28"D gives 



1 , 2i 



6(h)=2ih(T ) + — \h 2 -c h (34) 



and leads to a semicircle distribution for the h's 



p(h) = —^c h -h* (35) 
where the constant Ch is related to c by 

2 c . . 

c h = --—. (36) 

The strong coupling transition point occurs for the value of T at which p(h = 0) = 1, 
giving: 

c*=4=> T =p ~r = ^ == \ (37) 

c 



7T 



(tt — \JlX 2 — c) 



This can also be retrieved directly from equation(|32|). As a check we set c = so 
that M{9) is the unit matrix to give the partition function of QCD2 on the sphere 



and obtain the critical point discovered in [|nj which is 2T crit = ir , where 2T is the 
area. 

Another exactly solvable distribution a{6) = -j = ^ is studied in the appendix. 
This distribution is singular at the end points and provides an illustration of the case 
where there are several zeros to equation(|28|). 
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4 The Strong Coupling Phase 



In this section we present the analogue of the method described in JTT| for calculating 



the distribution of h's and the free energy in the strong coupling phase. In [11] the 
procedure was to start with the saddle point equation and look for the density p{h) 
which satisfies it subject to the constraint that for a finite interval [—b, b] p(h) takes 
on it's maximal value, ie. p(h) = 1 for he[— b, b\. The actual value of b is then found 
from consistency arguments. 

Although we derived our results from a completely different direction we still have 
a saddle point equation. It is the real part of equation (|30|) ie. 

P f dh'^-j = m(h)^ = 2Th-$Y, W (38) 

J II ll zeros 

The procedure is thus to solve the problem in the weak coupling phase by the inversion 
described above, and obtain if (/i) woak as an explicit expression of h, take the real 
part of this equation and solve (BH) for p(h) subject to the constraint that p(h) < 1. 



To do this we follow exactly the method discussed in |TTJ, which we repeat here 
almost verbatim. For simplicity's sake, we assume here that the distribution of /i's is 
symmetric and that there is only one interval, [—6, b], for which p(h) = 1. 

p{h) | e l sew here. 
Substituting into the real part of (53) we get 

3?F(/i) weak - log ^ = P J dh^j. (40) 

We define the resolvent for the distribution p(h) 

where the integral runs along the support of p(h), [—a, —b] and [b, a] (with a and b 
to be found later from the large h behaviour of H(h)), and generate the full analytic 
function H(h) from 3?iJ(/i) wcak by performing the contour integral 

H{h) = -J- J(a 2 -h 2 )(b 2 -h 2 ) £ ds ^ ( f )w ° ak = l ° 9 ~ b (42) 



2m 



(s - h)yj(a? - h 2 )(b 2 - h 2 ) 



where the contour encircles the cut of the square root but does not enclose any 
other cuts or singularities of the integrand. By inflating the contour we catch the 
pole at h = s, the cut of the log and any singularities of the analytic continuation 
of ?R,H(h) waak . The full resolvent H(h) stIons for the strong coupling phase involves 
adding to H(h) the contribution from the constant part of the full density p(h). This 
precisely cancels out the log term generated by the pole at h = s. We thus arrive at 
the expression for the resolvent for the strong coupling phase 

1 



H(h) strong = Wf(h)^ - J(a 2 -h 2 )(b 2 -h 2 ) I ds- 

v J-b I 



(s - h)^/(a 2 -h 2 )(b 2 -h 2 ) 
+ ( contribution from singularities of 9?iif(/i) weak ). (43) 

The last two terms give plus or minus in times the full density p(h). As is discussed in 



II] the second term can be written entirely in terms of the complete elliptic integral 
of the third kind EEf^, rrl: 

2b h-a 2Vabl 



lb — a 

7T 



(a + h){b + h) 



a + b h + 6 ' a + b 



(44) 



b + a\ (a - h)(b - h)' 

The coefficients a, and b are found by ensuring that for large \h\ H(h) stIong = \+0(-^). 
The large h limit of the elliptic integral is discussed in fLl] . 

For the semicircle distribution given as an example at the end of section (3) we 
have Kif(/i) wcak = —h. All the analysis conducted in |Ll| for the strong coupling phase 



is as before if we make the replacement in the notation of |TT| of A/2 — ► T/(l 



2T> 



5 Appendix 

5.1 An alternative approach 

Here we present an alternative method to finding the results of section(2). The 
charcter for U (N) is, up to a factor of i and a couple of trivial Vandermondt factors, 
the Itzykson Zuber determinant 

XR{e«) = AWJ^-i^^n^ ( 45 ) 
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where 1(0,^) = / dU H e Ntr ^ uH) -> e N2p ° M as iV -> oo. 



(46) 



Migdal in the context of induced QCD, and more recently Matytsin |14j] derived 
some useful results for the large N limit of the Itzykson Zuber integral: I[p(<p), cr{ip)} = 



,N 2 Fo[ P: 



CT 1 . We start from a useful formula which comes from taking the large N limit 



of |(^) 9 J(0,^). From Matytsin's paper |14| we quote the result (with the minor 
change of an extra factor of —i) 

1 1 



de'a(e')(-i9') q 



2-Kiq + l 



J dh' ( 



J + V(/i') + iirp(ti) 



dh 5p{h> 
/ 9 5F [p,cr 



+ V(/i') -mp(ti)) 9+1 



■ (4T) 



The integral runs along the support of p(h). We recognise V(/i) T ip{h) as the value 
of the resolvent H(h) above and below the cut and the equation simplifies to 

-1 1 



d6'(j{6'){6 



dh[iF(h) +iH(h)} 



9+1 



(4? 



27r q + 1 Jc 

where for convenience we have defined T = ^ ^jvy^ ; an< ^ ^ ne contour integration 
circles the cut. Multiplying both sides of the above equation by and summing 
from q = to oo we obtain 

dh 



e(6) = I ^ log[0 - iF(h) - iH(h)}. 
Jc 2n 



(49) 



This is the analogue of equation(^). It is an equation derived entirely from the 
properties of the character (Itzykson Zuber integral) itself. To proceed furthur one 
has to replace T{K) by its functional form derived from the saddle point equation 
evaluated for the particular action studied. It is important to remember that T is 
defined in terms of the Itzykson Zuber integral which is related to the character by 
the equation (|45|). The action has first to be written in terms of the Itzykson Zuber 
integral and then the saddle point equation evaluated. 

For the heat kernel we rewrite the action in terms of the Itzykson Zuber integral, 



dhp(h)A 2 (h)e~^^ h h(p,a) x A(i0)/(0) 



(50) 
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and then find the saddle point equation: 

T(h) = 2Th - 2V{h) (51) 

Substituting in to equation(pE9|) and using the identity H ± (h) — 2V(h) = —H T (h), 
where H T is the value of H(h) immediately above and below the cut, we obtain 

r dti 

0(9) = - 6 — log[0 - i2Th + iH(h)} (52) 
Jc 2n 

As earlier we take the derivative, expand the contour out to infinity, picking up the 
poles on the way, and arrive at the equation 

m m (53) 



zeros 



where h{9) satisfies 

9 -%2Th{9) + iH(h(0)) =0, (54) 



and the sum over zeros is the sum over all the solutions to equation (|54l) . We have 



thus found the analogue of equations (|30|) and(|28|). In the case where there are single 



solutions to equations (ESI) and (Bl), we see that equations Q53D and (Bl) are identical 



respectively, to equations ( p8|) and (|3lD - 

5.2 An example: The inverted semicircle distribution 

We present here the exact solution for a nontrivial density. The spectrum of eigen- 
values is 

9 k = gcos(—) (55) 
and the corresponding density and resolvent 

° {e) = 7W=V md °W = 7pW' (56) 



This density was used in the large N limit of the 2D chiral model studied in [15 . 



where it corresponded to an external magnetic field precisely tuned to excite the full 
spectrum of mass states of the theory. The heat kernel action investigated in this 
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article corresponds to the ID chiral model in an external magneitic field. Following 
the procedure of section(3), we solve equation(|28|). For convenience we introduce the 
rescaled variables; x = -^=, 77 = h\pXT and b 2 = in terms of which the equation 
is 

i 

71 + ix ; = (57) 

We square up the square root and arrive at the quartic 

x 4 - 2ir]x 3 - {rf + b 2 )x 2 + 2bir}X + {b 2 ri 2 - 1) = 0, (58) 



Although the quartic equation has 4 solutions only three satisfy the equation (|57|) . 
This can be intuitively understood by looking at r\ far from the origin. In this case 
there is one solution in which x is far from the origin (x ~ —iff) and the only two 
other possible solutions are the ones caught in the singularities of the end points of 
the square root. The final expression for the resolvent H(h) involves summing up all 
three of these solutions. Using a trivial identity x\ + £2 + ^3 + = 2b 2 ii] for the 
roots of a quartic allows us to express H(h) in terms of the single quartic solution 
(single sheet), x 4 , that does not satisfy the inversion equation (|57D. Using standard 
formulae for the roots of a quartic, the solution for the resolvent can be expressed as 

^) = -y|(^+*(v^+v^-v^)) ( 59 ) 



with zi, Z2 and z 3 defined by 



Zl 


= ^b 2 


- h 2 


+ l(F+ + F-) 


z 2 = 


¥- 


b 2 + 


1(ujF + + uj 2 F-) 


z 3 = 


¥- 


lv 2 + 


1(uj 2 F + + ujF~) 



(60) 
(61) 
(62) 



where lu — |(— 1 + iy/3) is the cube root of 1, and F ± are given by 



F ± = -\/k 3 - 18k + 54b 2 ± y (/c 3 - 18A; + 546 2 ) 2 - (k 2 - 12) 3 with k = r/ 2 + b 2 

(63) 
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All the cut structure resides in the square root term inside the definition of F . The 
\fz terms do not contribute cuts since for z ~ 0, z oc rj 2 , and the square root undoes 
the square. Similarly the cube roots do not contribute any cuts since, close to their 
zeros, the arguments are proportional to k 3 and the root undoes the cube. The zeros 
of the square root in equation(|63|) and thus the cutpoints for the quartic are 

V2= ±^-b 2 + ^ + (uB+ + u 2 B-) 

where uj — ~(— 1 + iV%) and B ± are defined by 

B± = ~^-jp + j + 8(3± ]p {p2 - 1)l with 13 = Vt 62 - (67) 

On the physical sheet (the sheet of the quartic corresponding to the resolvent) there 
is the single pair of cut points ±773. The other cut points connect together, in pairs, 
the three other non-physical sheets. The cut structure for H(h) is therefore given by 
a single cut or single interval for the support of the density running from —773 to 773. 

For b 4 < 7^ the support of the density lies entirely along the real axis, but for 
b 4 > ^ the end points move out into the complex plane. This is a signal that there 
is no longer a single Young tableau contributing to the sum. Probably it can be 
interpreted as the manifestation of the fact that a "density" of terms of alternating 
sign contributes in the large N limit rather than a single term. 

There are two cases where everything simplifies (the simplification reduces the 
inversion of equation([57]) to the solving of a quadratic). The first is when b = 0, 
which corresponds to all the eigenvalues lying at the origin, ie. when M(6) is the 
unit matrix. In this case both the cube root and the square root around the 2, can 
be taken explicitly and a semicircle distribution recovered as before. The other more 
useful case is when h=0. Again the cube roots and square roots around the z\ can be 
taken explicitly and we get a simple expression for the discontinuity across the cut 
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(64) 
(65) 
(66) 



at the origin ie 



^ + 4T 2 -f. (69) 
Since the density for the h's reaches it's maximum at the origin we have found the 



critical point for the distribution ie. p(h) max = 1 at T = f V 71 " 2 ~~ 9 2 - This result 
could also have been obtained directly from equation (|32|). 

Note 

Since the completion of this work it has come to our attention that D. Gross and 
A. Matytsin [16] have also been studying the large N limit of equation(Q) and have 



found, amongst other things, a similar simple inversion relation for the density of 
heighest weights for the case where g(0) = g{T)^ . 

Acknowledgements 

We would like to thank Andrei Matytsin for the simple expression for the strong 
coupling transition point. We are also grateful to Matthias Staudacher, Edouard 
Brezin and Michael Douglas for many fruitful discussions. 

References 

[1] J.S. Dowker, J. Phys. A (1970) 451. 

[2] P. Menotti and E. Onofri, Nucl. Phys. B190 (1981) 288-300. 

[3] A. A. Migdal, Sov. Phys. JETP 42 (1975) 413. 

[4] B. Rusakov, Mod. Phys. Lett. A5 (1990) 693. 

[5] D. Gross, Nucl. Phys. B400 (1993) 161. 

[6] D. Gross & W. Taylor, Nucl. Phys. B400 (1993) 181. 

[7] V.A. Kazakov & I.K. Kostov, Nucl. Phys. B176 (1980) 199. 

[8] V.A. Kazakov, Nucl. Phys. B179 (1981) 283. 

15 



[9] V.A. Kazakov & I.K. Kostov, Phys. Lett. 105B (1981) 453. 

[10] I.K. Kostov, Nucl Phys. B415 (1994) 29-70. 

[11] M.R. Douglas & V.A. Kazakov, Phys. Lett. 319B (1993) 219. 

[12] A. Matytsin, Private communication. 

[13] A.A. Migdal, Mod. Phys. Lett. A Vol8, No.4 (1993) p359-371. 

[14] A. Matytsin, Nucl. Phys. B411 (1994) 805. 

[15] V.A. Fateev, V.A. Kazakov & PB. Wiegman, LPTENS-94-07, hep- 
th@xxx.lanl.gov - 9403099 (1994). 

[16] D. Gross & A. Matytsin, PUPT 1503, hep-th@xxx.lanl.gov -9410059 
(1994) 



16 



